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The structures and properties of rhenium nitrides are studied with density function based first principle 
method. New candidate ground states or high-pressure phases at Re:N ratios of 3 : 2, 1:3, and 1 : 4 are 
identified via a series of evolutionary structure searches. We find that the 3D polyhedral stacking with strong 
covalent N-N and Re-N bonding could stabilize Re nitrides to form nitrogen rich phases, meanwhile, 
remarkably improve the mechanical performance than that of sub-nitrides, as ResN, Re2N, and Re3N2. By 
evaluating the trends of the crystal configuration, electronic structure, elastic properties, and hardness as a 
function of the N concentration, we proves that the N content is the key factor affecting the metallicity and 
hardness of Re nitrides. 



The combination of electron-rich 4d and 5d transition metals (TM) with light covalent-bond forming ele- 
ments (LE) like B, C, N, and O has been proved to be an effective way to design novel superhard material over 
the past decade^ l The idea is that the collaboration of TM's high electron concentration (EC) and the 
directional bonding formed by strong hybridizations between TM d electrons and LE s,p electrons can effectively 
withstand both elastic and plastic deformations. With this philosophy, compounds such as OsB2^, ReB2^, PtN2^'^, 
IrN2^'^, OsN2^, and PdN2^ were successfully designed. All of them are promising candidates for superhard and 
fracture-resistant materials with the elastic modulus comparable to diamond and c-BN. Recently, two Re sub- 
nitrides ResN (SG P-6m2) and Re2N (SG P6s/mmc) were synthesized by Friedrich et al}^ at moderate pressures 
(13-20 GPa) and temperatures (1600-2400 K). The ultra-high bulk modulus of > 400 GPa made them novel 
ultra-incompressible materials^^"^^. However, independent theoretical works by Zhang et al}^ and Zang et aV^ 
proved that the relatively lower ideal tensile strengths inherently limit the application of ResN and Re2N as 
superhard material. This discrepancy tells that the small compressibility is not equal to high hardness^°, and the 
high EC of Re sub -nitrides which is responsible for their low compressibility should actually go against the 
directional bonding. The Re:N stoichiometry may be the key to understand the competition between the EC and 
the directional bonding. And recent experimental and theoretical attempts to design N-rich Re nitrides are the 
subjects of intensive interest. Zhao et al?^ and Wang et al?^ proposed that Re mononitride forms in wurtzite-type, 
NbO-type, and NiAs-type structures, respectively. Besides, orthorhombic (SG Phcn) and tetragonal (SG P4/ 
mmm) Re dinitrides (denoted as orth-ReN2 and tetr-ReN2) with high bulk and shear modulus has been predicted 
theoretically^^"^^. Recently, Kawamura et aU^ synthesized the MoS2-type ReN2 with metathesis reaction. 
However, the MoS2-type ReN2 with a weak bulk modulus is less stable according to the calculations of Wang 
et a/.^^, and a monoclinic structure (SG C2/m) was proposed. In addition, the hypothetical ReP4-type ReN4 with 
high computational hardness has also been predicted by Zhao et aV^ and Aydin et al}^. 

Despite the substantial amount of works focused on Re nitrides. There are still two questions to be answered. 
Firstly, the multiple oxidation states of Re promise the emergence of Re-N phases at wide stoichiometrics and 
pressures. Therefore a comprehensive study to update the Re-N phase diagram toward different Re/N ratios and 
pressures, especially in the N-rich side, is necessary. Secondly, the effect of the N concentration on the mechanical 
behavior of Re nitrides and relevant mechanism are still unclear. To address these issues, we employ a series of 
evolutionary structure searches to reconstruct the low-T Re-N phase diagram at moderate pressures of 0- 
100 GPa. New stable phases for the known Re:N ratios and new stoichiometrics were updated and established. 
The trends of the crystal structure, electronic transition, and mechanical performance of Re nitrides as a function 
of N concentration are compared based on the updated full-scale Re-N phase diagram. The results show that the 
metallicity and hardness of Re nitrides can be correlated to the N content, and N-rich Re nitrides are better 
candidates for superhard materials due to the strong 3D covalent structural networks. 
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Table 1 | Space group (SG); calculated lattice parameters, a, b, and c (A); unit cell volume Vq (A^); and formation energy Ef(e V/atom) of 
Re, ResN, Re2N, Re3N2, ReN2, and ReNa at 0 GPa. And the parameters of high-P phase ReN4 at 50 GPa 

Phase SG a b c Vq Ef 



Ke 


1 his work 


Pds/mmc 


0 ~7~7 




4.4o 


1 A no 

1 4.yo 






D £ 0 0 

Rer. Jo 




i./O 




A A I. 

4.46 






1x63 IN 


This work 


r omz 


z .00 




7 1 A 
/ . 1 0 




kJ.kJO 




Ref . 1 0 




2.78 




7.15 








Ref. 1 8 




2.83 




7.19 






ResN 


This work 


Pds/mmc 


2.86 




9.87 


34.89 


0.03 




Ref. 1 0 




2.83 




9.88 








Ref. 1 8 




2.86 




9.88 






Re3N2 


This work 


P-6m2 


2.84 




8.30 


57.94 


-0.03 


ReN2 


This work 


C2/m 


6.83 


2.84 


9.39 


27.77 


0.19 




Ref. 27 




6.82 


2.84 


9.36 






ReNs 


This work 


Imm2 


5.25 


2.81 


4.75 


34.97 


0.44 


ReN4 


This work 


Cmmm 


7.54 


3.59 


2.63 


35.50 


0.83 



Results and discussion 

The global stability of Re nitrides can be quantified by constructing 
the thermodynamic convex hull within considered pressures, which 
is defined as the average atom formation enthalpy of the most stable 
phases at each composition: AH = H(ReJ<\y) - xH{Re) - yH(N2)/2. 
Any structure with the enthalpy lies on the convex hull is considered 
to be thermodynamically stable and synthesizable experiment- 
ally^°-^2. To give a more realistic prediction, both the experimental^^ 
and computationaP^'^^ high-P (low-T) phase orders of nitrogen were 
adopted (the C-N2 was replaced by 8-N2 in the calculations due to its 
structural uncertainty). The formation energy Ef of predicted Re 



nitrides can be calculated in the same way, listed in Table 1. The 
high-pressure convex hull and pressure- composition phase diagram 
of Re nitrides are depicted in Figure 1 and Figure 2. We present a rich 
phase content of Re nitrides toward different Re:N ratios and pres- 
sures. Particularly, three stable phases denoted as P-6m2-Re3N2, 
Jmm2-ReN3, and Cmmm-ReN4, beside the well-known RcsN, 
Re2N sub-nitrides, and N-rich ReN2, are proposed here as the com- 
plements to the Re-N phase diagram. It follows that the formation of 
Re nitrides, especially in the N-rich side, will be endothermic pro- 
cesses at absolute 0 K and 0 GPa due to their positive formation 
energies toward the thermodynamic summation of Re and a-N2. 
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0.8 1.0 0.0 0.2 0.4 

Concentration of N, x 

Figure 1 | Convex hull diagram for Re-N system. At pressures of (a) 10 GPa, (b) 40 GPa, (c) 70 GPa, and (d) 100 GPa, respectively. Solid circles 
represent stable phases while open circles denote metastable phases. N-rich stoichiometries are highlighted by cyan shadows. 
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Figure 2 | Pressure-composition phase diagram for Re-N system. 

Calculated with the (a) experimental (Ref. 33) and (b) computational 
(Ref. 34,35) phase orders of nitrogen, respectively. Red lines represent new 
phases proposed in this work. N-rich stoichiometrics are highlighted by 
cyan shadows. 

For Re sub -nitrides, ResN is stable at whole considered pressures 
while Re2N is thermodynamically stable only at high pressures, 
which is accord with the experimental results of Friedrich et al}^ 
and make our thermodynamic calculations reliable. In addition to 
the experimental sub -nitrides, our structure searches establish the 
ground state of new Re3N2 stoichiometry for the first time (denoted 
as P-6m2-Re3N2), which is thermodynamically stable up to 40.0 GPa 
as shown in Figure 1 and Figure 2. For the N-rich side, structure 
searching confirms that the C2/m-ReN2 proposed by Wang et al?'^ is 
the most stable phase and can be synthesized at the pressures of 
above 7.4 GPa. The previous orth-ReN2, tetr-ReN2, and MoS2-type 
ReN2^^ become metastable phases during the considered pressures. 
With further increasing the N concentration, a new stable structure 
(denoted as Imml-Ro^^) with extended Re:N ratio of 1 : 3 can be 
synthesized at wide pressures of 38.3-100 GPa (see Figure 2 (a)). 
And for Re tetranitrides, the proposed Cmmm-ReN4 is stable at 
elevated pressures of above 68.2 GPa, which can be viewed as the 
limit of high N content phase of the predicted 0-100 GPa Re-N 
phase diagram. Note that the hypothetical ReP4-ReN4 predicted by 
Zhao et al}^ and Ay din et al}^ should also be a metastable phase and 
can hardly be synthesized (see Figure 1 (c)). The adoption of the 
computational high-P phase order of nitrogen^^'^^ gives consistent 
phase diagram of Re nitrides, displayed in Figure 2 (b), except that 
the stable pressures of Jmm2-ReN3 and Cmmm-ReN4 have been 
raised to 39.5 and 89.8 GPa, respectively. 

The crystal structures of Re nitrides are schematically showed in 
Figure 3. The mechanical and dynamical stabilities are examined by 
calculating the elastic constants (see Supplementary Table SI online) 
and phonon dispersions (see Supplementary Figure SI online), all 
proposed phases are mechanically stable at 0 GPa with the elastic 



constants meet the Born- Huang stability criteria^^. Besides, the P- 
6m2-Re3N2 is dynamically stable at 0 GPa. And N-rich C2/m-ReN2 
and Jmm2-ReN3 are tested to be dynamically stable from 0 GPa to 
their synthetic pressures, which means that they are quenchable to 
the ambient condition. The Cmmm-ReN4 is dynamically stable at 
high pressures (50-100 GPa). According to Figure 3, a unique struc- 
tural transition induced by increasing N content can be clearly 
revealed. The Re-Re-N-Re, Re-N-Re, and Re-Re-N-N-Re lamellar 
stacking orders of Re3N, Re2N, and Re3N2 sub -nitrides has been 
broken and transformed into intriguing polyhedral stacking config- 
urations as in C2/m-ReN2, Imm2-ReN3, and Cmmm-ReN4. Among 
them, the structure of Re3N2 is special since it can be viewed as the 
substitution of the isolated N with N2 dimer in the configuration of 
Re3N. This reveals that the segregation of Re layers by N bilayer can 
stabilize Re nitride to low energy and also weak metallicity (see the 
following electronic structure discussions). With this idea, an ana- 
logous structure can be built for ReN that is by substituting the N 
with N2 dimer in Re2N along the c-axes. Along with slightly stagger 
of the neighbouring two Re layers, it results in P-3m-ReN, which is 
more stable than the wurtzite-ReN, NbO-ReN, and NiAs-ReN pro- 
posed by Zhao et al?^ and Wang et aU^ (see Figure 1 (a)), and can 
easily be obtained by the 1 : 1 fixed Re:N ratio evolutionary structure 
searches. However, the lamellar order with a growing number of the 
vertical N2 dimers destabilize Re nitride. The P-3m-ReN is meta- 
stable and will decompose into Re3N2 and ReN2 as shown in 
Figure 1. Besides, the introduction of N2 dimer can lower the crystal 
density and results in large volume phases. The PV works will sac- 
rifice the enthalpies markedly, which is responsible for the instability 
of Re3N2 at high pressures. The polyhedral stacking configuration 
can be an effective way to stable Re nitride and forms in dense N-rich 
phases. For C2/m-ReN2, the stacking consists of the corner or edge 
shared ReNy decahedrons and also bonded with N-N connections 
according to Figure 3 (e). With an additional N into the lattice, the 
stacking for Jmm2-ReN3 consists of more regular ReNy octahedrons 
with the connections of the special tri-nitrogen units (see Figure 3 
(f)). As to the high-pressure phase Cmmm-ReN4, the novel ReNg 
cubic box can actually be formed due to the coordination number 
of Re changing from 7 to 8. The parallel ~ 1 .40 A N-N bonds (0 GPa) 
in C2/m-ReN2, Imm2-ReN3, and Cmmm-ReN4 represent the typical 
single N-N bond lengths in PM nitrides, e.g., 1.41 A for PtN2, 1.43 A 
for OSN2, and 1.42 A for IrN2, therefore we can draw the conclusion 
that N-N forms in single bond since the bond-length similarity 
within a related bonding scheme^^. The common polyhedral stacking 
configurations of N-rich Re nitrides is a prerequisite for superhard 
since it can avoid the direct Re- Re interactions, and the structural 
networks can be constructed primarily by Re-N and N-N covalent 
bonds. The space group, calculated equilibrium lattice parameters, 
and the unit cell volumes of rhenium nitrides with various stoichio- 
metries are listed in Table 1. Our calculated lattice parameters of Re, 
Re3N, Re2N, and C2/m-ReN2 are all in good accordance with pre- 
vious studies^^'^^'^^'^^ 

The updated full-scale Re-N phase diagram allows us to evaluate 
the trend of the mechanical performance of Re nitrides as the func- 
tion of N concentration, and correlate it with the structural and 
electronic transitions. The calculated bulk modulus (5), shear modu- 
lus (G), Young's modulus (£), Poisson's ratios (v), densities (p), and 
hardness (Hy) of Re nitrides at different N concentrations are 
depicted in Figure 4 and Table 2 (for comparability, the elastic data 
of the high-P phase Cmmm-ReN4 in Figure 4 are achieved at 0 GPa), 
while the electronic density of states (DOS) are showed in Figure 5. 
Re nitrides have similar metallic bonding features owing to the finite 
DOS at the Fermi level (£p), which originate mostly from 5d electrons 
of Re. According to Figure 4 (a) and (e), the high N content in the 
crystal can decline the bulk modulus and densities markedly, which 
heralds the decrease of the incompressibility and EC of N-rich Re 
nitrides. This is accords with the reducing of the metallicity, i.e., the 
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Figure 3 | Crystal structures of Re nitrides, (a) Re; (b) Re^N; (c) RcsN; (d) P-6m2-Re3N2: Rel atoms occupy the 1/(0.6667, 0.3333, 0.5), Re2 atoms 
occupy the 2h (0.3333, 0.6667, 0.2361), andN atoms occupy the 2i (0.6667, 0.3333, 0.0841) positions; (e) C2/m-ReN2; (f) Imml-ReN^: Re atoms occupy 
the 2b (0.5, 0, 0.3466), Nl atoms occupy the 4c (0.21, 0.5, 0.2455), and N2 atoms occupy the 2b (0, 0.5, 0.4123) positions; and (g) Cmmm-ReN^ (at 
50 GPa): Re atoms occupy the 2b (0, 0.5, 0), while N atoms occupy the 8^ (0.8401, 0.1927, 0.5) positions. Large and small balls denote Re and N atoms, 
respectively. 



DOS at the Fermi level 0<l(Ep)), of Re nitrides at elevated N concen- 
trations. Here, we should note that the Imml-ReN^ with the lowest B 
possess the N(£p) minimum of —0.1 1 eV/atom (see Figure 5). On 
the other hand, a parallel increasement of shear modulus, Young's 
modulus, and a reduction of Poisson's ratios can be clearly revealed 
with increasing N contents, as shown in Figure 5 (b), (c), and (d). The 
unsmooth change of the G, £, and v between the Re^N, Re2N, and 
Re3N2 sub -nitrides can be attributed to the discontinuity of the C12 
and Ci3 elastic constants (see Supplementary Table SI online), which 
relate to the Poisson effect in the hexagonal structures. According to 
Table 2, the calculated mean value of the G and E for N-rich C2/m- 
ReN2, Imm2-ReN3, and Cmmm-ReN4 are about 46 and 86 GPa 
respectively higher than that of Re and sub-nitrides. In addition, 
the Poisson's ratios of Jmm2-ReN3 and Cmmm-ReN4 have been 
dropped down from the typically delocalized metals (—0.3) to the 
level of strong directional covalent materials (—0.2). The reinforce- 
ment of the N-rich Re nitrides has its own origin and can be attrib- 
uted to the enhanced hybridizations between N 2p and Re 5d 
electronic states (i.e., in the bonding states of —5 to 0 eV below the 
Fermi level) as shown in Figure 5. Therefore, we can conclude that 
the primary contribution of the high N content to N-rich Re nitrides 
is that it facilitates more p electrons involved in bonding, thus can 
translate the material from high EC strong metal sub -nitrides to the 
strong covalent N-rich nitrides with enhanced directional bondings. 
These results suggest that the N content is a key factor to influence 
the hardness of Re nitrides and the hardness of C2/m-ReN2, Imm2- 
ReN3, and Cmmm-ReN4 should be far higher than that of Re3N, 
Re2N, and Re3N2 sub-nitrides. 

The hardness of Re nitrides as a function of N concentration is 
estimated by recently proposed empirical modeP^ with the formula 
Hy = 2.0(k^Gy '^'-3.0, as shown in Figure 4 (f). It can be seen that the 
high N content promote the hardness of Re nitrides according with 
the changes of G, £, and v as discussed above. On the other hand, the 
strong covalency of the bonding in N-rich Re nitrides allows us to 
evaluate their intrinsic hardness by microscopic hardness mod- 
els^°"^^ Here we use the method proposed by Gao et al^° and its 
further development which includes a small metallic component of 
chemical bond and considering the orbital form of s-p or s-p-d^^. The 



hardness of the ^-type particular bond in metallic and s-p-d hybri- 
dized crystals is written as^°'^^'^^: 

H(f(GPa) = AoiN^fidT'-'e- ''si/," -32.2(f:)« " ( 

where Aq is the proportionality coefficient takes 350 and 1051 for the 
pure covalent N-N and polar covalent Re-N bonds, respectively. is 
the valence electron density of the bond, d is the bond length. The 
Phillips iconicity/^ of chemical bond is calculated by the Mulliken 
population^^. The/^ in Equation (1) is a factor introduced to account 
for the metallicity effect on hardness and can be calculated as = 
kBTDp/rie, where tie is the total number of valence electrons in the 
unit cell, Dp is the density of electronic states at the Fermi leveP^. 
The total hardness can be expressed as the average of all bond hard- 
nesses: „ 

H,= n(HO"" (2) 

where n^' is the number of bond of type /n composing the actual 
complex crystal. The detailed bond parameters of hardness calcula- 
tions can be found as Supplementary Table S2 online. The hardness 
of C2/m-ReN2 and Jmm2-ReN3 at ambient condition are evaluated 
as 38.4 and 44.4 GPa respectively, which could be potential candi- 
dates to be superhard materials. 

Conclusion 

In summary, we have systematically explored the structures and 
properties of Re nitrides at wide stoichiometries and pressures with 
ab initio evolutionary method. We corrected the stable structure of 
ReN4 to the phase with the Cmmm symmetry, and new stoichiome- 
tries of P-6m2-Re3N2 and Jmm2-ReN3 are established for possible 
synthesis. It is found that high-pressure can be an effective approach 
to stable Re nitrides. The C2/m-ReN2, Imm2-ReN3, and Cmmm- 
ReN4 can be synthesized at pressures above 7.4, 38.3, and 68.2 
GPa, respectively. And C2/m-ReN2 and Jmm2-ReN3 are retainable 
to ambient condition. A good relation between the mechanical beha- 
vior and N concentration of Re nitrides was established by detailed 
evaluating the variations of the crystal configurations, electronic 
structures, elastic properties, and computational hardness. Our 
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Figure 4 | Elastic properties of Re nitrides at 0 GPa. Calculated (a) bulk modulus B (GPa), (b) shear modulus G (GPa), (c) Young's modulus £ (GPa), 
(d) Poisson's ratio v, (e) density p (g/cm^), and (f) Vickers hardness Hy (GPa) of Re nitrides as a function of N concentration. The neighboring data points 
were connected by B-Spline to evaluate the changing trends. N-rich stoichiometries are highlighted by cyan shadows. 



results show that the high N content benefits the formation of the 
covalent 3D polyhedral stacking structure and can effectively avoid 
the weak Re-Re metal interactions, which result in the reinforcement 
of the Re nitrides from ultra-incompressible sub -nitrides to super- 
hard N-rich materials. These results are expected to stimulate the 
exploration and discovery of new material, and the approach of 



harden Re nitride by high LE content may also applicable to other 
TM compounds. 

Methods 

Evolutionary structure searches are performed with USPEX code"^^"*^. Structure cells 
with Re:N ratios of3: 1,2:1, 1.5:1, 1 : 1, 1 : 1.5, 1:2, 1:2.5, 1 : 3, 1 : 4, 1 : 5, and 1 : 6 



Table 2 


Calculated bulk modulus 


8 (GPa), shear 


modulus G (GPa 


, Young's 


modulus 


E (GPa), Poisson's ratio v, density p (g/ 


cm^), and 


Vickers hardness Hy (GPa) of Re-N system at different pressures 












Phase 




SG 


P 


B 


G 


E 


V 


P 


Hv 


Re 


This work 


Pds/mmc 


OGPa 


376 


185 


477 


0.29 


20.72 


15.49 


ResN 


This work 


P-6m2 


OGPa 


397 


202 


518 


0.28 


19.12 


17.26 




Ref . 1 0 






395 














Ref. 1 8 






397 


203 










ResN 


This work 


P62/inmc 


OGPa 


399 


189 


490 


0.30 


18.40 


14.98 




Ref. 1 0 






401 














Ref. 1 8 






401 


182 










Re3N2 


This work 


P-6m2 


OGPa 


379 


209 


529 


0.27 


16.82 


19.63 


ReN2 


This work 


C2/m 


OGPa 


376 


210 


531 


0.26 


12.81 


19.86 




Ref. 27 






369 


217 










ReNs 


This work 


Imm2 


OGPa 


330 


234 


568 


0.21 


10.84 


29.51 


ReN4 


This work 


Cmmm 


OGPa 


354 


283 


670 


0.18 


10.15 


38.72 


ReN4 


This work 


Cmmm 


50 GPa 


528 


379 


919 


0.21 


1 1.33 


40.82 



SCIENTIFIC REPORT! | 4 : 4797 j DOI: 1 0.1 038/srep04797 



5 



TDOS 



<Re-d 



N-/7 



.Re(OGPa) ^ 
\ i , 


1 ' 1 . 


-7^-6w2-Re.N(0GPa) 


— i ' 
1 
1 

; ■ ^r— 


-P6//7?/7?6-Re N (0 GPa) 


: A 



g 4 _ /^-6w2-Re (0 GPa) n Ni i 



g 2 - C7/m-ReN^ (0 GPa) 
Q 

4 




2 |-/wm2-ReN3 (0 GPa) 
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within 2 and 4 formula units (f.u.) are implemented at 0, 30, 60, and 100 GPa, 
respectively. To update a full-scale Re-N phase diagram, moreover, variable-com- 
position structure predictions^" are also performed at 0, 50, and 100 GPa, respect- 
ively. The first generation of structures was produced randomly and the succeeding 
generations were obtained by applying heredity (60%), atom transmutation (10%), 
and lattice mutation (30%) operations. Each structure was fully relaxed to an energy 
minimum at different pressures with the Vienna ah initio simulation package 
(VASP)^\ The projector augmented waves (PAW) method^^ within the generalized 
gradient approximation (GGA) of Perdew-Burke-Ernzerhof (PBE)^^ was adopted. By 
performing accurate convergence tests, we choose a cutoff of 520 eV for the wave 
function expansion into plane waves and a k-mesh of 0.03 X 27i A"^ within the 
Monkhorst-Pack scheme for sample the Brillouin zone, which ensures the error bars 
of total energies are less than 1 meV/atom. Besides, all forces on atoms were con- 
verged to less than 1 meV/A and the total stress tensor was reduced to the order of 

0. 01 GPa. Phonon curves were calculated by the direct supercell method with the 
forces obtained by the Hellmann-Feynman theorem^"^'^^. Elastic constants were cal- 
culated by the strain-stress method, while bulk modulus B, shear modulus G, Young's 
modulus £, and Poisson's ratio v are derived from the Voigt-Reuss-Hill averaging 
scheme^*^. The theoretical Vickers hardness was estimated by using the Chen's 
modeP'^ and the Gao's modeP°, respectively. For the latter, the Mulliken bond 
population is calculated with the CASTEP code^^'^*^. The ultra-soft Vanderbilt pseu- 
dopotential (USPP)^^ with PBE-GGA is chosen. The cutoff energy and k-mesh are 
same as the setting in VASP calculations. 
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